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Abstract 

We perform a state-of-the-art global fit to all Higgs data. We synthe- 
sise them into a 'universal' form, which allows to easily test any desired 
model. We apply the proposed methodology to extract from data the 
Higgs branching ratios, production cross sections, couplings and to 
analyse composite Higgs models, models with extra Higgs doublets, 
supersymmetry, extra particles in the loops, anomalous top couplings, 
invisible Higgs decay into Dark Matter. Best fit regions lie around the 
Standard Model predictions and are well approximated by our 'uni- 
versal' fit. Latest data exclude the dilaton as an alternative to the 
Higgs, and disfavour fits with negative Yukawa couplings. We derive 
for the first time the SM Higgs boson mass from the measured rates, 
rather than from the peak positions, obtaining Mh = 124.2 ± 1.8 GeV. 



1 Introduction 



After the discovery of a new particle around 125.5 GeV announced by the ATLAS [1] and 
CMS [2] LHC collaborations during 2012, all LHC and TeVatron collaborations presented at 
the Moriond 2013 conference the new results based on the full collected data. These include 
the most important 77, ZZ* and WW* channels as well as updates to the fermionic channels. 
Such results will stay with us for next two years until LHC with full energy starts operating. 
Therefore it is the right moment to analyse their implications. 

We want to know if the new particle is the long-waited Standard Model (SM) Higgs boson 
[3, 4, 5, 6]. On one side, the experimental collaborations are measuring its discrete quantum 
numbers to check if it is a scalar. On the other side, various theoretical groups [7, 8, 9] started 
to approximatively reconstruct from data its production cross section and its decay modes 
and consequently its couplings to check if they agree with the SM predictions or with other 
models beyond the SM. Clearly, this is a more significant test that can be precisely done only 
by the experimental collaborations, which indeed started to present analyses along these lines. 
However these experimental fits, presented in the form of likelihood plots within a few specific 
beyond-the-SM models, are of little use to theorists who are interested in different models. 

We here propose how experimental collaborations could report their results in a model- 
independent and useful way, such that these results would be readily and reliably used by 
theorists who want to test any desired model. The new ingredient that we introduce and that 
allows for this simplification is the assumption that new physics can be approximated as a 
first-order perturbation with respect to the SM predictions. We find that this assumption is 
increasingly supported by data, that agree with the SM with precisions around the 20% level. 

Such results, obtained after two years of LHC operation and with only 25/fb data per 
experiment, implies severe constraints on models where the Higgs boson is a portal to new 
physics. We analyse several models and rule out alternative scenarios to the Higgs boson. 

The paper is organised as follows. In section 2 we present the data and our fitting procedure. 
In section 3 we derive the first measurement of the Higgs mass from the rates, rather than from 
the position of the peaks in the 77 and ZZ invariant mass distributions. In section 4 we 
present the 'universal' format for data mentioned above. Next, in section 5 we present fits 
in various specific models, updating our previous results and comparing the full fit to the 
simplified 'universal' fit to verify that it is a good approximation. We fit Higgs cross sections in 
section 5.1, Higgs couplings in 5.2, composite Higgs models in 5.3, new physics in loops in 5.4, 
two Higgs doublet models in 5.5, the MSSM in 5.6, the dilaton in 5.7, the Higgs invisible width 
in 5.8 and models where DM couples to the Higgs in 5.9. In section 6 we summarise the results 
and draw our conclusions. 

2 The data 

Searches for the SM Higgs boson have been carried out in proton-proton collisions at yfs = 7 
(2011 data) and 8 TeV (2012 data) with about 25/fb of total integrated luminosity. 
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There are four main production modes for Higgs boson from pp collisions at ^/s ~ 8 TeV. 
The gluon-gluon fusion production mode has the largest cross section, followed in turn by vector 
boson fusion (VBF), associated Wh and Zh production, and production in association with top 
quarks, tih. The cross sections for the Higgs boson production modes and the decay branching 
fractions, together with their uncertainties, are taken from [10]. 

Our updated analysis uses the new data presented at the Moriond 2013 conference by the 
CMS, ATLAS and TeVatron collaborations [11, 12, 13, 14] in the following five decay modes: 
77 [15], ZZ* (followed by ZZ* decays to U, 2£2u, 2£2q, 2£2t) [16], WW* (followed by WW* 
decays to tvtv^tvqq) [17, 12], t + t~ (followed by leptonic and hadronic decays of the r-leptons) 
[13, 18] and bb [13, 19], and the first tentative measurements in the /x 4 ' \i~ [20] and Z7 channels 
[21], as well as their combination [22]. Here and throughout, i stands for electrons or muons 
and q for quarks. 

For a given Higgs boson mass, the search sensitivity depends on the production cross section 
of the Higgs boson, its decay branching fraction into the chosen final state, the signal selection 
efficiency, the mass resolution, and the level of standard model backgrounds in the same or a 
similar final state. For low values of the Higgs boson mass, the h — >• 77 and h — > ZZ* — > A£ 
channels play a special role due to the excellent mass resolution for the reconstructed diphoton 
and four-lepton final states, respectively. The h — > WW* — > Ivtv channel provides high 
sensitivity but has relatively poor mass resolution due to the presence of neutrinos in the final 
state. The sensitivity in the bb and t + t~ decay modes is reduced due to the large backgrounds 
and poor mass resolutions. 

We include in our data-set all exclusive 77 and rr sub-categories described by the experi- 
mental collaborations by telling how much each Higgs production channel in the SM contributes 
to the various rates. We adopt the Multivariate Analysis (MVA) 77 analysis from CMS and we 
combine all experiments, such that we find an average 77 rate very close to the SM prediction. 
Consequently our results differ from previous analyses performed without including the latest 
CMS 77 data [7]. 

This is an important issue because, while most of the presented LHC results are well con- 
sistent with the SM predictions within experimental errors, there are few unexpected new 
developments that need commenting. The most important of them is the discrepancy between 
the ATLAS and CMS results in the h — > 77 channels. With full luminosity, ATLAS finds an 
overall rate of 1.65 ± 0.34, higher than the SM prediction of 1, and higher than the CMS result 
of 0.80 ± 0.30. The two measurements are compatible within 2a. In addition, the two CMS 
77 analyses (MVA and cut based) show different signal rates. Finally, the two Higgs boson 
mass determinations in ATLAS, from the peaks in the 77 and ZZ channels, differ by 2a. Both 
experiments have cross checked their analyses and reached conclusions that those deviations 
are due to statistical fluctuations of both signal and background. This conclusion implies that: 
(i) combining all data in a global fit is meaningful and increases the precision; (ii) selecting 
instead any single measurement, for example the ATLAS excess in 77, is not justified and 
introduces a bias in the data. 

The experimental collaborations report Higgs boson rates R in units of the central value 
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Figure 1: Measured Higgs boson rates at ATLAS, CMS, CDF, DO and their average (horizontal 
gray band at ±la). Here (red line) corresponds to no Higgs boson, 1 (green line) to the SM 
Higgs boson (including the latest data point, which describes the invisible Higgs rate). 

of the SM prediction. Their results could be fully encoded in a likelihood C(R,Mh), but only 
a limited amount of information is reported by the experiments. Sometimes the experimental 
collaborations report the measured rates as R exp ± R e ". Sometimes collaborations only report 
the upper bounds on rates at 95% C.L., -Ro™e rve d; an d the expected upper bound at 95% C.L. in 
absence of a Higgs boson signal, -Rejected; as function of the Higgs boson mass m h . Assuming 
that the y 2 = —2 In £ has a Gaussian form in R, these two experimental informations allow 
one to extract the mean R exp and the standard deviation R eri as i? exp = -^observed ~ -^expected 
and i? err = -Rejected/ 1-96, where 1.96 arises because 95% confidence level corresponds to about 
2 standard deviations. The x 2 is approximated as 



y 2 = V - 1)2 (1) 



where the sum runs over all measured Higgs boson rates I. 

The theoretical uncertainties on the Higgs production cross sections <jj start to be non- 
negligible and affect the observed rates in a correlated way. We take into account such cor- 
relations in the following way. We subtract from the total uncertainty Rf r the theoretical 
component due to the uncertainty in the production cross sections, obtaining the purely exper- 
imental uncertainty, R e / T ~ exp . The theoretical error is reinserted by defining a x 2 which depends 
on the production cross sections, 
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and marginalising it with respect to the free parameters (jj, constrained to have a central value 
cr* h and an uncertainty <T° rr given by 

<r(pp ->• h) th = (19.4 ± 2.8) pb, c(pp -> jj/i) th = (1.55 ± 0.04) pb, 

a{pp ->■ W7i)th = (0.68 ±0.03) pb, <r(pp -> Z/i) th = (0.39 ± 0.02) pb, . (3) 

a{pp ->■ tt/i) th = (0.128 ± 0.018) pb at ^/s = 8TeV. 

We summarise all data in fig. 1 together with their la error-bars. The grey band shows the 
±lcr range for the average of all rates: 1.00±0.10. It lies along the SM prediction of 1 (horizontal 
green line) and is about 10 a away from (the horizontal red line is the background-only rate 
expected in the absence of a Higgs boson). 



3 Reconstructing the Higgs mass 

The CMS and ATLAS collaborations reported measurements of the pole Higgs mass Mh 
obtained as the position of the peaks observed in the invariant mass of the h — > 77 and 
h — > ZZ — > 4£ distributions: 



M h = 125.66 ± 0.34 GeV 



125.4 ± 0.5 stat ± 0.6 syst GeV CMS 77 

125.8 ± 0.5 stat ± 0.2 syst GeV CMS ZZ 

126.8 ± 0.2 stat ± 0.7 syst GeV ATLAS 77 

124.3 ± 0.6 stat ± 0.5 syst GeV ATLAS ZZ 



(4) 



These measurements are mutually compatible, and the uncertainty is so small that in the 
subsequent fits to rates we can fix M h to its combined best-fit value. 

Indeed, we here discuss how the Higgs mass can be independently measured, with a bigger 
uncertainty, by requiring that the measured rates agree with their SM predictions. Such pre- 
dictions have a dependence on the Higgs mass that, around 125 GeV, can be approximated as 



a(pp -> X) w a(pp -> X) Mh=125GeY x [1 + c x x (M h - 125 GeV)]. (5) 

In table 1 we list the values of the coefficients cx and of the measured rates for the various 
processes averaging all experiments, as well as the Higgs mass indirectly derived from such 
rates. We see that the single best indirect determination of comes from the h — > WW 
rates, that presently have no sensitivity to Mh if one wants to measure it from a mass peak. 
On the other hand, the h — > 77 signal that offers the best peak measurement of Mh has very 
little indirect sensitivity to Mh, because the 77 rate happens to have a weak dependence on 
Mh- Averaging over all channels we find 

M h = 124.2 ±1.8 GeV (Higgs mass extracted from the rates, assuming the SM) (6) 

which is compatible with the determination of the pole Higgs mass obtained in a model- 
independent way from the positions of the peaks. 
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Process X 


h -> WW 


ZZ 


h — > 77 


Vh ->■ 


/l — ^ TT 


Sensitivity cx 
Measured rate/SM 
Higgs mass in GeV 


6.4%/ GeV 
0.84 ±0.17 
123.0 ±3.0 


7.8%/ GeV 
1.06 ±0.22 
126.2 ±2.7 


-1.5%/ GeV 
1.07 ±0.19 
121 ± 12 


-5.4%/ GeV 
1.19 ±0.42 
122 ±8 


— 4.1%/ GeV 
1.11 ±0.28 
123 ± 7 



Table 1: Determinations of the Higgs mass from the rates, assuming the SM. 

4 The universal Higgs fit 

We perform the most generic fit in terms of a particle h with couplings to pairs of t, b, r, W, Z, g, 7 
equal to r t , r b , ry, r^, r ff , r 7 in units of the SM Higgs coupling. This means, for example, that 
the coupling to the top is given by r t (m t /V)hit, where r t = 1 in the SM and V = 246 GeV is 
the Higgs vacuum expectation value. Similarly, the /177 coupling is assumed to be ry times its 
SM prediction. In the SM this couplings first arises at one loop level. Experiments are starting 
to probe also the hfin and the hZ'j effective couplings, so that also the corresponding and 
rz-y parameters will start to be measured. This discussion can be summarized by the following 
effective Lagrangian: 



171+ , - rrih , - m T , m T , M% , 9 2Mw , 

U = r t ^htt + r b ^hbb + r T ^hfT + r^hw + r z ^hZl + r w —^hW+W- + 

^77 a u tp w t_ „ „99 as un a n a 1 ». ^ z 7 " 



The various SM loop coefficients csm are summarised in appendix A. This Lagrangian is often 
written in a less intuitive but practically equivalent form by either using SU(2) L £g> U(l)y- 
invariant effective operators, or assuming that the Higgs is the pseudo-Goldstone boson of a 
spontaneously broken global symmetry and writing its chiral effective theory [7]. 

Furthermore, we take into account the possibility of Higgs decays into invisible particles X 
(such as Dark Matter or neutrinos) with branching ratio BR inv . In almost all measured rates 
BR inv is equivalent to a common reduction r of all the other Higgs couplings, BR inv ~ 1 — r 2 , 
such that BR inv is indirectly probed by data [8]. The only observable that directly probes an 
invisible Higgs width is the pp —> Zh — > XX rate measured by ATLAS [30] , which implies 

BR inv = -0.19 ±0.43. (8) 

Any possible new-physics model can be described as specific values of the r» parameters. 
Several examples are provided in section 5. 

Following the procedure described in the previous section, we approximatively extract from 
data the function 

X 2 (r t , r b , ry, r w , r z , r g ,r^, r Zl , ?>, BR inv ) (9) 

which describes all the information contained in Higgs data. We find \ 2 = 58.8 at the best fit 
(56 data points, 10 free parameters), marginally better than the SM fit, x|m = 61-7 (no free 
parameters). 
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Figure 2: x 2 as function of the model-independent Higgs couplings to the various SM particles, 
varying them one-by-one. 

4.1 Universal fit to small new physics effects 

The universal x 2 °f ec L- (9) nas a too complicated form to be reported analytically, and depends 
on too many variables to be reported in numerical form, like plots or tables. For these reasons, 
previous analyses [7, 8, 9] focused on particular BSM models with a reduced number of param- 
eters. For example, fig. 2 shows the fit as function of each r i; setting all others to their SM 
values of unity: we see that the x 2 are approximately parabolic. 

We here observe that Higgs data are converging towards the SM predictions with small 
errors, thereby it is time to start making the approximation 

Ti = 1 + Cj with 6i < 1 (10) 

and BR inv = e inv . The observable rates Rj are computed at first order in e iy and consequently 
the x 2 is expanded up to second order in e^. As well known, this Gaussian approximation 
is a great simplification; for example marginalizations over nuisance parameters just becomes 
minimisation, which preserves the Gaussian form. Fig. 2 suggests that this approximation 
already seems reasonably good. 

For LHC at 8 TeV the main observables are approximated as 

Rh-+ww = 1 - 1.14e 6 + 1.58e 9 

Rhr+zz = l-1.14e 6 + 1.58e fl 

R h ^ TT = 1 - 1.14e 6 + 1.58e ff 

i4-> 77 = 1 - 1.14e6 + 1.83e ff 

R h ^ bb = 1 + 0.86e 6 + 1.58e 9 

Rv(h^bb) = l + 0.86e fe -0.17e 3 

The full x 2 can now De reported in a simple form. Indeed the x 2 is a quadratic function of the 
€i, and it is usually written as 

X 2 = J^(e; - ^i){(J 2 )^^j - Mi), where (a 2 )^ = on'.jO , (12) 

i,3 



- e inv - 0.17 e t + 1.72e w + 0.02e z 

- e inv - 0.17 e t - 0.28evK + 2.02e z 

- e inv + 1.83 et - 0.28e w + 0.02e z 

- e inv - 0.18 et - 0.456^/ - 0.06e z + 2e 

- e inv - 0.17e t - 0.28e w + 0.02e z 

- e inv - 0A8e t + 0.83e w + 0.67e z . 
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in terms of the mean values fii of each parameter of its error Oi and in terms of the correlation 
matrix py. We believe that this is the most useful form in which experimental collaborations 
could report their results. From our approximated analysis of LHC and TeVatron data we 
obtain: 
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-0.05 ±0.34 
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-0.09 ±0.15 


P = 




0.48 
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0.43 


0.20 


0.31 
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±0.03 ±0.19 
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0.48 


0.34 


0.36 


0.61 


0.54 
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We have not reported the central value of r t = 1 ± e t , of e^ 7 and of e M because they presently 
are known only up to uncertainties much larger than 1. Future searches for tih production, for 
h — > and for h — > fi + /i~ will improve the situation. 

In many models the Higgs couplings to vectors satisfy ew = £z, because of SU(2) l invariance. 
Furthermore, in many models LEP precision data force ew and ez to be very close to 0. 
This restriction can of course be implemented by just setting these parameters to be equal or 
vanishing in the quadratic y 2 . 

Since the uncertainties on the ej parameters are now smaller then 1, the universal approxi- 
mation starts to be accurate. In the next sections, where we analyze several specific models, we 
will systematically compare our full numerical fit (plotting best fit regions in yellow with con- 
tinuous contours at the 90 and 99% C.L.) with the universal approximation (best fit ellipsoidal 
regions in gray with dotted contours, at the same confidence levels). 



5 Model-dependent Higgs fits 

5.1 Higgs production cross sections 

Assuming the SM predictions for the Higgs decays, we extract from the data the Higgs pro- 
duction cross sections. We find that they agree with SM predictions, as shown in fig. 3a. As 
expected, the most precisely probed cross section is the dominant one, a(pp — y h). At the op- 
posite extremum a(pp — > jjh) is still largely unknown. The uncertainties on the reconstructed 
cross sections are correlated, altought we do not report the correlation matrix. 

5.2 Higgs couplings 

We here extract from data the Higgs boson couplings to vectors and fermions, assuming that 
only the SM particles contribute to the h — > gg, 77, 7Z loops. This amounts to restrict the 
universal fit in terms of the parameters by setting the parameters for loop couplings to 

r g = r t , r 7 » 1.282^ - 0.282^ r Zl » 1.057r w - 0.057r t (14) 
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Figure 3: Left: reconstruction of the Higgs production cross sections in units of the SM pre- 
diction. Right: reconstruction of the Higgs couplings to the t, Z,W,b,r, assuming that no 
new particles exist. The SM predicts that Higgs couplings are proportional to particle masses 
(diagonal line). 

These numerical expressions are obtained by rescaling the expressions for the SM loops sum- 
marised in appendix A. In particular, the W loop (rescaled by r w ) and the top loop (rescaled 
by r t ) contribute to h — > 77 with a negative interference. 

Under this assumption the top coupling of the Higgs, r t , becomes indirectly probed via the 
loop effects. The fit to the couplings is shown in fig. 3a and agrees with the SM predictions, 
signalling that the new boson really is the Higgs. We do not report the correlation matrix, 
which can be immediately obtained by inserting eq. (14) into the universal x 2 °f ec L- (12). 

5.3 Composite Higgs models 

Models where the Higgs is composite often assume the further restriction, in addition to eq. (14), 
of a common rescaling with respect to their SM values of the Higgs boson couplings to the W, Z 
bosons and a common rescaling of the Higgs boson couplings to all fermions. These rescalings 
are usually denoted as a and c, respectively: 

r t = r b = r T = r^ = c, r w = r z = a. (15) 

The resulting fit is shown in fig. 4a. We see that the best fit converged towards the SM; in 
particular data now disfavour the solution with c < which appeared in previous fits and 
gave an enhanced h — > 77. We also see that our approximated universal fit (dotted contours) 
reproduces very well our full fit (continuous contours). 
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Figure 4: Left: fit of the Higgs boson couplings assuming common rescaling factors a and 
c with respect to the SM prediction for couplings to vector bosons and fermions, respectively. 
The two sets of contour lines are our full fit (continuous) and our approximated 'universal' fit 
(dotted). Right: values of the y 2 along the trajectories in the (a, c) plane shown in the left 
panel, and given by a = y/1 — £ and c = a (magenta) c = (1 — 2£)/a (blue) c = (1 — 3£)/a 
(red), as motivated by composite Higgs models. The black dashed curve corresponds to a = 1 
and c — 1 — £. 



In fig. 4b we show the full \ 2 restricted along the trajectories in the (a, c) plane (plotted in 
the left panel) predicted by simple composite pseudo-Goldstone Higgs models in terms of the 
parameter £ = (V/F n ) 2 , where is the scale of global symmetry breaking. 



5.4 New physics only in the loop processes 

We here assume that only the loop processes are modified with respect to the SM predictions, 
summarized in appendix A. This amounts to restrict our universal fit setting 



r t = H = r T = r^ = r w = r z = 1, r — = r r — = r 16 

r(/i ^ ^)sm 5 r(/i -»■ 77)sm 



with BR inv = and r Zl = 1. The latter assumption is at present justified because of the large 
experimental error in the h — > rate, even though in general new physics in the loop processes 
would induce deviation from unity in both r^ 7 and r 7 . The result is shown in the left panel of 
fig. 5, under the form of a fit to the ratios of BR(/i — > gg) and BR(/i — > 77) with respect to the 
SM. One can see that the SM is well within the la contour. The universal fit approximates well 
the full fit. The dashed trajectories show the loop effect due to extra scalar particles with the 
same quantum numbers of the top (red), of the bottom (blue), of the tau (vertical black line). 
The explicit expressions for the contribution of scalar, fermion and vector particles running 
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0.0 o.5 1.0 1.5 2.0 0.0 0.5 1-0 1.5 2.0 2.5 3.0 

BR(/i->gg)/SM Scalar coupling to the Higgs r s 



Figure 5: Left: fit for the Higgs boson branching fraction to photons and gluons, with 1 and 2a 
contours. The dashed curves shows the possible effect of extra scalar partners of the top (red), 
of the bottom (blue), of the tau (black). Dotted lines show the Gaussian approximation. Right: 
Upper bound at 90% (solid) and 99% (dashed) C.L. on the new scalar coupling rs to the Higgs 
as a function of the new scalar mass ms- 

in the loop can be found in appendix A. Note that any additional color-less but electrically 
charged particle would lead to the same trajectory obtained for the scalar partner of the r. 

To better investigate the constraints on a possible new scalar 5*, in the right panel of fig. 5 
we show the upper bound, as function of the scalar mass ms, on the scalar coupling r s to the 
Higgs boson, defined by the coupling 

rs^hSS. (17) 

The resulting loop effects are summarised in appendix A. The solid and dashed curves in fig. 5b 
are respectively the upper bounds at 90% (solid) and 99% (dashed) C.L. More stringent limits 
are obtained on the top and bottom partners than on the r partner. 

One can also use the universal fit with the assumption of eq. (16) to derive indirect con- 
straints on the top quark magnetic and chromomagnetic dipole moments [26, 27], which in the 
SM are expected to be respectively g t ~ 2 and k t ~ 2. Allowing g t and k t to vary freely, the 
h — > 77 and h —> gg amplitudes are modified with respect to the SM according to: 



( 

where the quantities and are defined in eq. (33) of the Appendix. Fig. 6 shows the 



r — — — V8 1 2/ r - -k 2 - - (-\R) 

7 r (w) , (t) ' 9 8 ' 2 ' { } 
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Top (chromo)magnetic dipoles 




90,99% CL 



1.6 1.8 2.0 2.2 2.4 

Top chromomagnetic dipole k, 



Figure 6: Best fit regions for the (chromo) magnetic dipole moments of the top quark. 

90% and 99% CL. allowed regions for g t and k t . The uncertainty on kt is comparable to the 
one from its direct measurement, while the one for g t is even smaller [28]. 

5.5 Models with two Higgs doublets 

There are four types of two Higgs doublets models (2HDM) where tree level flavour-changing 
neutral currents (FCNCs) are forbidden by a Z 2 symmetry [35] and both doublets Hi and H 2 
get a vacuum expectation value: 

• type I [36, 37] where only one doublet couples to all quarks and leptons; 

• type II [37, 38], where up-type quarks couple to H 2 and Hi couples to down-type quarks 
and leptons. The Higgs sector of the MSSM is a type II 2HDM; 

• type X (lepton-specific or leptophilic) where H 2 couples only to quarks and Hi couples 
only to leptons; 

• type Y (flipped) [41], where H 2 couples to up-type quarks and H 2 to down-type quarks, 
and (contrary to the type II HDM) leptons couple to H 2 . 

For an extensive review see [39] and for some previous fits see [40]. The modification to Yukawa 
couplings to up-type and down-type quarks and leptons in the four 2HDMs are: 

Type I Type II Type X (lepton-specific) Type Y (flipped) 

r t cos a/ sin (3 cos a/ sin (3 cos a/ sin (3 cos a/ sin (3 

r b cos a/ sin (3 — sin a/ cos (3 cos a/ sin (3 — sin a/ cos f3 

r T cos a/ sin (3 — sin a/ cos (3 — sin a/ cos (3 cos a/ sin (3 
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Figure 7: Fit to the t-quark and to b-quark and r-lepton Yukawa couplings assuming the struc- 
ture predicted by the various types of two Higgs doublet models. The point marked as C SM' is 
the Standard Model; the point marked as 'FP ? is the fermiophobic case. 
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As usual, tan/3 = V2/V1 is the ratio of the VEVs of the two doublets and a is the mixing angle 
of the CP-even mass eigenstates. The SM limit corresponds to — a = tt/2. In all of the 
models the vector couplings are also modified as 

r w = r z = sin(/3 — a). (19) 

The results of our fits are presented in fig. 7 in terms of the fermion couplings r t , r b , r T , restricted 
by the 2HDM models to lie within the green regions. 

The type II 2HDM (upper panel) allows for independent modification of the t coupling r t , 
and for a common modification of the b and r couplings, r b = r T . The former is predicted 
be reduced and the latter enhanced by the model. The modification of eq. (19) of the vector 
couplings can be equivalently written as r w = r z = (l+r t r b )/(r t +r b ) ~ l + e t e b /2, showing that 
it is a small second order effect. In this model a negative t Yukawa coupling is still allowed at 
slightly more than 99% CL. The red line in the same panel shows the parameter space allowed 
by type I 2HDM, where all the couplings scale uniformly. 

In the flipped 2HDM (middle panel) the r Yukawa coupling changes in the same way as 
the t coupling and the region with negative coupling is disfavoured by data. Finally, in the 
leptophilic 2HDM (lower panel) the t and b couplings vary in the same way, while the r coupling 
is independent. 

The universal fit provides a good approximation to the full fit in all 2HD models. 



5.6 Supersymmetry 

Supersymmetry can affect Higgs physics in many different ways, such that it is difficult to make 
general statements. We here focus on the two most plausible effects: 

• Stop loops can enhance or reduce the Higgs coupling to the top (and consequently the 
h gg, 77, Zj rates) by an amount given by 



Rt 



m. 



■ ■ v + 

1 + f 



2 



+ 



ti 



2 

m- 
12 



(A t - /i/tan/3)^ 

2 2 

tl t2 



(20) 



in the limit of heavy stop masses, mi 3> m t . 



• The type II 2HDM structure of supersymmetric models modifies at tree level the Higgs 
couplings, as already discussed in section 5.5. 

All of this amounts to specialise the universal x 2 inserting the following values of its parameters 



cos a 



sin a 



r w = r z = sin(/3 — a) 



(21) 



r t = R t^-15i r b = r T = r^ = -, 

sin p cos p 

Furthermore, the parameters r 9 ,r 7 ,r^ 7 relative to loop processes are fixed as in eq. (14). We 
trade the a parameter (mass mixing between Higgses) for the pseudo-scalar Higgs mass m^ 
using 

m\ + Mf 



tan 2a 



m\-M\ 



tan 2/3. 



(22) 
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Figure 8: Left: Fit to the two main effects present in supersymmetry: stop loop correction to 
the hti coupling and tree-level modification of the Higgs couplings due to the two-Higgs doublet 
structure. Right: fit as function of the (3 -function coefficients 63 = 6 7 that parameterise dilaton 
models. The SM Higgs is reproduced at the experimentally favored point b 3 = 6 7 = ; while the 
pure dilaton is excluded at more than 5a. 

Finally, we assume a large tan j3, as motivated by the observed value of the Higgs mass. Fig. 8a 
shows the resulting fit. Once again, the universal fit approximates well the full fit. Of course, 
supersymmetry can manifest in extra ways not considered here, e.g. very light staus or charginos 
could enhance h —> 77 [31]. 

5.7 Data prefer the Higgs to the dilaton 

As another example of a model where both the tree level and the loop level Higgs couplings 
are modified, we consider the dilaton. The dilaton is an hypothetical particle (p, that, like the 
Higgs, couples to SM particles with strength proportional to their masses [32]. More precisely 
the dilaton has a coupling to the trace of the energy-momentum tensor T^, suppressed by some 
unknown scale A: 




The dilaton couplings to gg and 77 differ from the corresponding Higgs boson couplings, because 
eq. (23) contains the latter two quantum terms, that are present in T£ because scale invariance 
is anomalous and broken at quantum level by the running of the couplings. Indeed 63 and 6 7 
are the /3-function coefficients of the strong and electromagnetic gauge couplings. In the SM 
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0.0 0.2 0.4 0.6 0.8 1.0 10 20 30 40 50 60 70 



Invisible Higgs BR DM mass in GeV 

Figure 9: Left: fits to the invisible Higgs boson branching fraction under the two different 
assumptions described in section 5.8. The full fit (continuos curves) is well approximated by 
the universal fit (dotted curves). Right: upper limit on the spin-independent DM cross section 
on nucleons as a function of the DM mass for scalar (green), Majorana fermion (red) and 
vector (blue) DM. We adopted the 95% C.L. bounds BR inv < 0.19 (solid, eq. (25)) and < 0.28 
(dot-dashed, eq. (26)). The shaded region is excluded at 90% C.L. by XenonIOO [25]. 

they have the explicit values 63 = —7 and fe 7 = 11/3: we call 'pure dilaton' this special model, 
which gives a significant enhancement of h -B- gg. 

Models where a dilaton arises usually often contain also new light particles, such that 63 
and b 1 can differ from their SM values. Thereby we perform a generic fit where 63 and b 1 are 
free parameters in addition to A. Then, our universal fit is adapted to the case of the generic 
dilaton by setting 

V 

r = r w = r z = r t = r b = r T = — , r g « r(l — 1.45fe 3 ), r 7 « r(l + 0.156 7 ) (24) 
where V = 246 GeV. 

In our previous analyses [8, 9], the dilaton gave fits of comparable quality to the SM Higgs, 
despite the significantly different predictions of the dilaton: enhanced 77 rates and reduced 
vector boson fusion rates. The first feature is no longer favoured by data, and the second 
feature is now disfavoured: so we find that present data prefer the Higgs to the 'pure dilaton' 
at about 5a level. We then consider the generic dilaton, showing in fig. 8b that the allowed part 
of its parameters space is the one where it mimics the Higgs, possibly up to a sign difference 
in r g and/or r 7 . The dilaton becomes identical to the SM Higgs in the limit 63 = 6 7 = 
and A = V. This situation is not easily realisable in models, given that adding extra charged 
particles increases 6 7 rather than reducing it. 
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5.8 Higgs boson invisible width 

Next, we allow for a Higgs boson invisible width, for example into Dark Matter. We perform 
two fits. 

1. In the first fit, the invisible Higgs width is the only new physics. We find (blue curves in 
fig. 9a) that present data imply 

BR inv = -0.08 ± 0.16 i.e. BR inv < 0.19 at 95% C.L., (25) 

2. In addition to the invisible width we also allow for non-standard values of h — > 77 and 
h — > gg, finding a weaker constraint on BR inv (red curves in fig. 9a) 

BR inv < 0.28 at 95% C.L. (26) 

Notice that the main constraint con BR inv does not come from the direct search for pp — > Zh — > 
U$ T (included in our data-set) but from the global fit [8, 33]. 



5.9 Dark Matter models 

The invisible Higgs boson decay width [33] constrains Dark Matter (DM) candidates with mass 
below Mh/2. The Higgs sector of the SM allows for a direct coupling to particles of a hidden 
sector. If the latter are stable and interact weakly with the SM sector, they could represent 
viable Dark Matter (DM) candidates. If DM particles have mass below Mh/2, the Higgs boson 
can thus decay into a couple of DM particles, which would escape detection. Invisible Higgs 
decays are constrained by the fact that the ATLAS and CMS Higgs rates are compatible with 
the predictions of the SM Higgs boson. The experimental bound on BR inv can be used to 
constrain the DM mass and its elastic cross section on nucleons probed in direct detection 
experiments, as illustrated for instance in [24], where DM is assumed to be either a scalar S, 
or a Majorana fermion / or a vector V coupled to the Higgs as 

rs^hSS + r f ^fhff + r v 2 -^hV,V, . (27) 

The partial Higgs decay width into dark matter T(h — > DM DM) and the spin-independent 
DM-proton elastic cross section ergi can be calculated in terms of the parameters of the above 
Lagrangian. Both are proportional to the square of the DM-Higgs coupling, so that the ratio 
fi = cr SI /r(/?, DM DM) depends only on the the unknown DM mass and on the known masses 
and couplings of the relevant SM particles (see for instance the expressions provided in [24]). 

This allows us to relate the invisible Higgs branching fraction to the DM direct detection 
cross section: 

T(h — »■ DM DM) _ ctsi 
rS M + r(/i DM DM) ~ /iTl M + a S i 

where = 4.1 GeV is the total Higgs decay width into all SM particles, that we fix to its SM 
prediction. For a given DM mass, an upper bound on the Higgs invisible branching fraction 
implies an upper bound on the DM scattering cross section on nucleons. The relation between 



RR _ 1 V 1 ~Z ^ 1V1 ^ L ^> - " (oa\ 
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the invisible branching fraction and the direct detection cross section strongly depends on the 
spinorial nature of the DM particle, in particular, the strongest (weakest) bound is derived in 
the vectorial (scalar) case. 

Imposing the upper bounds on BR inv derived in section 5.8, fig. 9 shows the corresponding 
upper limits on the spin-independent DM cross section on nucleons as a function of the DM 
mass, in the case of scalar (green), Majorana fermion (red) and vector (blue) DM candidates. 

In all cases, the derived bounds are stronger than the direct one from XENON100 as long 
as the mass of DM is lighter than This conclusion does not rely on the assumption that 

DM is a thermal relic that reproduces the observed cosmological DM abundance. 

6 Discussion and Conclusions 

The LHC experiments reported their measurements of Higgs boson properties at the Moriond 
2013 conferences, based on the full collected luminosity during 2011 and 2012. At the same 
time, TeVatron reported their final Higgs results. With the crucial inclusion of the full CMS 
77 data (missing in previous analyses), at this stage all main Higgs results from TeVatron 
and from the first phase of LHC have been basically presented. Those results will drive our 
understanding of particle physics, until new 13 TeV LHC data will be available. 

Motivated by these results, we have performed a state-of-the-art global fit to Higgs boson 
data, including all sub-categories studied by the experimental collaborations, for a total of 56 
data-points, as summarised in fig. 1. We found that the average Higgs rate is 1.00 ± 0.10 
in SM units, supporting the SM Higgs boson hypothesis. The Higgs boson mass is usually 
determined from the peaks in the invariant mass distribution of ZZ and 77. We performed the 
first measurement of the Higgs boson mass from the rates, finding that the two determinations 
are compatible: 



The LHC physics program has been successful: with only « 25/fb of data per experiment 
the Higgs boson has been discovered and several of its properties determined within ±20% 
precision. We are now entering into the era of precision Higgs physics — deviations from the SM 
due to new physics no longer can dominate the data. This observation allowed us to propose a 
'universal' form in which experiments could report their results allowing theorists to easily test 
any desired model. The new assumption that makes possible this significant simplification is 
that new physics is a small correction to the SM. While we present our own global combination 
in 'universal' form in eq. (13), we stress that only the experimental collaborations can perform 
a fully precise analysis. 

We studied several new physics scenarios beyond the SM. We determined from data the 
production cross sections (assuming standard Higgs decays) and the Higgs decays widths (as- 
suming standard productions), finding that they lie along the SM predictions. In a more general 
context, we allowed all possible Higgs boson couplings to any SM particle to deviate from its 
SM value, finding that couplings to the W, Z, t, b, r must lie around their SM predictions up to 




125.66 ± 0.34 GeV from the peaks, 
124.2 ± 1.8 GeV from the rates. 



(29) 
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uncertainties of about ±20% (see fig. 3b). In particular, non-standard Higgs boson couplings to 
vectors, predicted by composite Higgs models, are most stringently constrained. The scenario 
of negative Higgs coupling to fermions ('dysfermiophilia') that gave the best fit with early LHC 
data is now disfavoured at more than 2a. 

We considered various specific new physics models: new scalars, 2HDM, supersymmetry, 
dilaton, composite Higgs, invisible Higgs decays, possibly into Dark Matter particles, anomalous 
couplings of the top, etc. The results of those fits are presented in numerous figures throughout 
the paper. Qualitatively, all reach to the same conclusions: 

i) best fit regions lie along SM predictions, imposing constraints on new physics; 

ii) our simple 'universal' approximation to the full fit is adequate. 

In particular we find that, with the latest data, the dilaton alternative to the Higgs is now 
excluded at 5<r, with the exception of the special non-minimal dilaton tuned to exactly reproduce 
the Higgs (section 5.7). 

We will update this paper at the light of future results. 

Acknowledgement This work was supported by the ESF grants 8943, MTT8, MTT59, 
MTT60, MJD140 by the recurrent financing SF0690030s09 project and by the European Union 
through the European Regional Development Fund. The work of P. P. G. has been partially 
funded by the "Fondazione A. della Riccia". 

A New physics contributions to loop processes 

The coefficients in the second line of eq. (7) arise at one-loop. They are obtained by summing the 
contributions of all scalars (S) fermions (/) and vectors (V) that couple to the Higgs as in eq. (27). 
The explicit expressions for the loop effects are [29] : 

c f ] = ^rsAsirs) = 2C f 2 T f A f {r f ) (30) 

c (5) = ^ rsAs{Ts) c Cf) = ^lr f A f (y) cOO = -^TyAyM 

where for each particle p = S, f, V, t p = rn^/Anip, N p is the number of colors, Cf is the Casimir of 
the color representation (Tr (T a T b ) = C25 ab ), and the loop functions are 

Mr) = ^ if(r) ~ r] , A f (r) = l)/(r) + r] (31) 

A v (t) = ^ [3(2t - l)/(r) + 3r + 2r 2 ] (32) 

with /(t) = arcsin 2 (-y/r) for r < 1 such that A p (t p ) — > 1 in the limit t p — > (heavy p-particle). 

In particular, in the SM, the hgg coupling is dominated by the top loop, and the /177 coupling 
arise from the sum of the top and W boson loops: 

41 = cf = A f (r t ) 41 = C W + cf ) = 2 Q A f (r t ) - -A v (r w ) . (33) 
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Beyond the SM (BSM) physics affects the parameters r g and r 7 as 

r.99 



1 + r 7 = l + ta. (34) 



C SM C SM 



For example, additional scalar particles with the same quantum numbers of a stop, sbottom and stau 
respectively contribute to Cgg M and to Cgg M as: 



' ' "" V -,As(t,) <•;," = () 

(35) 



cf = ftAsW cf = - A r~ b A s ( n ) cP = 
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